Spontaneous chirality via long-range electrostatic forces 
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We consider a model for periodic patterns of charges constrained over a cylindrical surface. In 
particular we focus on patterns of chiral helices, achiral rings or vertical lamellae, with the constraint 
of global electroneutrality. We study the dependence of the patterns' size and pitch angle on the 
radius of the cylinder and salt concentration. We obtain a phase diagram by using numerical and 
analytic techniques. For pure Coulomb interactions, we find a ring phase for small radii and a chiral 
helical phase for large radii. At a critical salt concentration, the characteristic domain size diverges, 
resulting in macroscopic phase segregation of the components and restoring chiral symmetry. We 
discuss possible consequences and generalizations of our model. 
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The ubiquitous role electrostatics play in the organiza- 
tion of biomolecular assemblies is a main pillar of modern 
biophysics (see, for instance, [i] and references therein). 
The co-assembly of oppositely charged molecules of- 
ten occurs over cylindrical structures with directional de- 
pendent functionalities, such as in fibers or filamentous 
viruses [i 3| for signal transmission or cellular structure 
mobility |4| . Other important cationic-anionic cylindrical 
co-assemblies include actin-binding proteins structures 
5 1 ^cylindrical micelles Q , and peptide-amphiphile fibers 
71, [sl] . A considerable aspect of these systems is the pres- 
ence of surface charge heterogeneities or domains [9, 10], 
which arises from the competition of electrostatic forces 
with segregating forces such as steric, van der Waals, and 
hydrogen bonding interactions. These surface patterns 
have been shown to be important for the stability and 
functionality of the aggregates [l, 11, 12. [l3|. 

A key feature in many biological processes is pattern 
recognition and specificity within isomeric aggregates 
[l3 | . Biomolecules exploit repetitive patterning over sur- 
faces, which often breaks some underlying symmetry, in 
order to increase their functionality [15|. In this con- 
text, chiral symmetry (that is the absence of improper 
rotation axis) is notably interesting: although its role in 
chemistry and biology is widely recog nized SSEIll, 
its origin in biology is the subject of much debate 18 1. 
Many systems are chiral due to the chirality of its com- 
ponents 3, 13] ■ For instance, colloidal solutions of rod- 
like viruses provide an intriguing example of chiral struc- 
tures in liquid crystalline phases (2L] . Nevertheless, un- 
derstanding whether electrostatic interactions alone are 
capable of producing chiral systems would shed further 
light on the issue. 

From a theoretical vantage point, the electrostatic pat- 
terning of a system of charges on cylindrical surfaces is 
however relatively unexplored. Certainly, the effects of 
long-range electrostatic forces have been widely studied 
for planar two dimensional systems [2^; also the behav- 
ior of short-range interactions over cylindrical geometries 



has been addressed, such as the Ising model on the cylin- 
der [l^]. We analyze here an intermediate case where 
charges are confined over a cylindrical surface and inter- 
act via long-range forces. A further issue we consider is 
whether spherically symmetric electrostatic interactions 
are capable to break translational, rotational or chiral 
isometrics of the cylinder. Recently, there has been in- 
terest to study crystalline systems over constrained ge- 
ometries such as the surface of spheres, cylinders, and tori 
24, 25, [2^. The generalization to more general curved 



substrates shows an interesting rich behavior [27'] . 

In this Letter we provide the full phase diagram of 
lamellar charged patterns on a cylinder, as a function of 
the cylinder radius and screening length. We explicitly 
show the existence of a phase where the system sponta- 
neously adopts chiral configurations. We study also its 
"Coulombic" origin by observing how it disappears when 
increasing the screening length. Our model is a gener- 
alization of the one introduced by F.J. Soils et al. 
for charge domains on cylindrical cationic-anionic co- 
assemblies in aqueous solution. We consider a 1:1 sto- 
ichiometric mixture that covers the surface completely. 
The surface can generally be considered incompressible 
due to the strength of hydrophobic interactions. The 
degree of the incompatibility among the cationic and an- 
ionic components - characterized by the standard Flory- 
Huggins parameter x ~ drives the molecules to separate 
and minimize the interfacial energy between the positive 
and negative surface groups. While x leads the system 
to segregation, the long-range Coulomb potential favors 
mixing. Molecules of like charge form homogeneous do- 
mains of average charge density a; boundaries between 
such domains acquire a line tension 7. In the strongly 
segregated limit we assume the lowest-energy configu- 
ration corresponds to domains arranged on a periodic 
lattice {A}. This system can be described by a coarse- 
grained two- terms free energy T = Ti+ T2 of a Wigner- 
Seitz cell, with the line tension term J-i opposing the 




FIG. 1: Schematic view of an unwrapped cylinder with a 
generic hehcal pattern, showing the notation used in the text. 

surface charge density term J-2 
^1=7^, d^^d^y<j{^)U{^-y)<j{y), (1) 

where I is the interfacial length between the domains, 
(t(x) is the local charge density, and JJ' represents an in- 
tegration over the whole surface and the other restricted 
to a single Wigner-Seitz unit cell. In real systems the 
electrostatic potential is attenuated by free ions in the 
solvent and we adopt a screened Coulomb potential: 
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where e is the dielectric constant, and k is the inverse 
Debye-Hiickle length. The latter takes effectively into ac- 
count possible salt concentration in the solution, and we 
will use it here also as a tunable parameter to control the 
long-range/short-range nature of the interactions. The 
two terms in the free energy T generate a characteris- 
tic length Lq = ^/7e/(T, and it is therefore convenient to 
adopt dimensionless units. All lengths x are measured in 
units of Lq, X = x/Lq, I = I/Lq, the inverse screening 
length becomes k = kLq, the charge density a — cr/i§, 
and the free energy is measured in units oiJ-o = 7^0, i-e. 
J- — T/Tq. For the sake of compactness in the rest of this 
Letter we assume that the name of a variable stands for 
its dimensionless counterpart, and we drop all the tilde 
symbols. We consider the free energy density F = T j 
where is the area of the Wigner-Seitz cell. A self- 
consistent minimization of the free energy F for a planar 
system gives lamellar patterns at 1:1 stoichiometric ratio 
[9| , and it has been confirmed by direct numerical simula- 
tions [2^. Henceforth we concentrate on the behavior of 
such lamellar patterns and their critical properties on the 
cylinder. In that case, we have F\ = 2/L where L is the 
distance between two neighboring lamellae [§] . The dou- 
ble integral in Fi can be computed first by "unwrapping" 
the surface of the cylinder of radius R onto an infinite set 
of parallel stripes of size I-kR on the plane, and then by 
summing all the pairwise electrostatic interactions over 
the planar periodic lattice {A} of the lamellae, by Fourier 
decomposition method. To simplify the calculations we 
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FIG. 2: Phase diagram of our model. Region I is the ring 
phase; Region II is the helical phase; Region III is the axial 
phase; the region at the top is the fully segregated phase. The 
solid lines are discontinuous transitions and the dashed line 
is a continuous crossover. 



adopted a square lattice oriented along the lamellae (see 
Fig. [T|). Formally: 

^2 = ^E/ / rf'/'^(0'^(A + ry)F(A-fry-C), 
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with D — y 4i?2 sin^ {px/{2R)) + pi is the projected dis- 
tance of p = (pxiPz) on the {x, z} plane, and 9{x) is the 
step function. By introducing the reciprocal lattice {Q}, 
defined by Q • A = 27rm, m G Z, and by using the Poisson 
summation formula we have: 



^2 = ^E'^(QM-Q)^(Q)' 
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where (t(Q) and T^(Q) are the Fourier transforms of a 
and V, respectively. They are: 



V{Q) = Rlg^RiOKg^RiO^ 
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where <^ = R^ + q1 and Ii, and Ki, are modified Bessel 
functions of order v of the first kind and second kind, 
respectively. For simplicity we used the twofold nota- 
tion Q — {qx,<lz) in the reference system parallel to the 
lamellar domains, and Q = {qx,Qz) in the reference sys- 
tem parallel to the axis of the cylinder (see Fig. [T]). Let 
Q be the pitch angle between the axis of the cylinder z 
and the direction of the lamellar stripes z. Any lattice on 




FIG. 3: Evolution of the pitch angle Q of helical lamellae 
with the Debye screening inverse-length k, at different cylin- 
der radii R. 



the surface of a cylinder must be commensurate with its 
circumference, that is it must be periodic in the x direc- 
tion with period 2t:R. The commensurability constraint 
in the reciprocal space read — 2'Km/L or equivalently 
Qj. = 2'KmcosQ / L and = 27rrnsinO/L, m S Z. More- 
over qxR = mn must hold, where n is the number of 
lamellar stripes per pitch. 

We numerically minimized the free energy density F with 
respect to all possible lamellar lattices, i.e. with respect 
to the spacing L and the pitch angle Q for several val- 
ues of R and k,. Our results are collected in the phase 
diagram in Fig. [2l There are four regions: a phase I 
characterized by ring patterns, a chiral phase II charac- 
terized by helical patterns, a phase III characterized by 
vertical lamellar patterns, and a fully segregated phase. 

Phase I is limited to fibers with circumference smaller 
than the domain size. This is natural when consider- 
ing the limit R <^ L where cosQ ~ nL/i-wR 1, i.e. 
Q ~ 7r/2, which is precisely the observed ring phase. 

Phase II is the primary result of this work. It proves 
the existence of chiral configurations (both right-handed 
and left-handed) due to electrostatic long-range interac- 
tions. A surprising feature of this phase is that at k = 
there is a preferred pitch angle Q* ~ tt/5 (see Fig. [3]). 
The long-range nature of electrostatic interaction may 
be the reason of such a preference, since by increasing 
K > the pitch angle decreases continuously down to the 
vertical lamellar phase III. Our numerical findings can 
be further justified on a theoretical ground by looking at 
the large- i? asymptotic expansion of the free energy. In 
the expansion Qx is kept fixed while the density (t(x) is 
not dependent on R since it is defined locally over a unit 
cell. For the potential V^(Q) we obtain (up to 0(l/i?'^) 
terms): 

V(o) 1 (4Q^ ~ 5gg - (g,^ + .^) 1 

(8) 

where Q = \/ q^^ + 9^ • The leading term is the planar 
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FIG. 4: By increasing the screening length above a critical 
value Kc, the free energy global minimum jumps from a finite 
L to 1/ = 00. In this plot R = 1000 and = 1.37 ± 0.05 

limit studied in ^28]. It is straightforward to show that 
the minimum of the free energy at k = with re- 
spect to the pitch angle is where dF(Q)/d9 = 0. From 
eq. ([8]) at K = we find 6* = arccos \/3/5 ~ 0.68 rad 
which is close to the numerical value 7r/5 ~ 0.63 rad. 
For K > the first-order correction to O* is negative and 
scale as k} . This is in agreement with the numerical be- 
havior depicted in Fig. [31 that the pitch angle initially 
decreases when increasing k. From eq. ([8]) we obtain also 
the value kq — 2^j2ji-K j L at which the pitch angle van- 
ishes, i.e. where phase II merges continuously into III 
and the helical patterns evolve into vertical stripes. We 
find L = 5.8 ± 0.2 at large- i? (see Fig. 2]), and therefore 
kq = 0.9 ± 0.2. The nature of phase II can further be 
seen in the discontinuous transition to the achiral phase 
I. Assuming the separation distance L is not strongly 
dependent on i?, there is a critical ratio LfR such that, 
as R decreases, discontinuously jumps from O* to 7r/2 
(see Fig.©. 

All the three phases we considered so far are bound at 
large k by a fully-segregated phase. Such a salt-induced 
phase transition is a first-order phase transition as the 
free energy jumps from a finite cell size to an infinite cell 
size (see Fig. |4]). Let Kc{R) be the critical value of k at 
which this transition occurs. In Fig.[3]the transition is at 
the jump of the pitch to 7r/2, near the "melting" transi- 
tion point K ~ Kc- At large- i? we find Kc = 1.6±0.1. The 
melting transition can be interpreted by a scaling argu- 
ment. At large-L the free-energy scales as F w 
If K is small then the global minimum is at a finite value 
of L, that is the size of the Wigner cell (see continuous 
black line in Fig. [4]). If k is large then the minimum 
jumps to infinity, that is to a macroscopically segregated 
system (see dotted green line in Fig.H]). 

We note that the phase diagram for the cylinder is 
richer than its planar counterpart. In fact, due to the 
commensurability constraints of the cylindrical surface 
new lamellar pattern configurations with different chiral 
angles arise for different radii R and lamellar spacing L. 
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The theoretical prediction for the cyhndrical geometry in 
the large radius limit agrees well with the result for pla- 
nar geometry. However a striking distinction from the 
planar case is the appearance of a preferred chiral angle 
of about &* = arccos a/3/5 ^ tt/S, which is numerically 
found to be conserved over a wide range of radii. It is an 
amusing coincidence this angle corresponds to the one 
found in recent studies of capsid-protein arrangements 
on cylindrical viruses J[2l| as well as in biomimetic self- 
organized nanotubes [2]. However Coulomb interaction 
is not in principle the driving force of the assemblies con- 
sidered in 0, [ilj. Therefore an interesting question is 
whether such a pitch angle is universal, and what is its 
origin. We attempt an answer by considering the large- i? 
asymptotic expansion for higher-order decaying isotropic 
interactions, V3 ^ l/D^. We find (up to 0{l/R^) terms: 
^(Q) ~ —T^Q + TTql/iSQ^R"^), whose stationary points 
are only at 8 = and tt/2 (this statement holds true 
also for higher multipolar terms). Although this argu- 
ment cannot be considered as proof, it shows that the 
universality of the pitch angle is unlikely. However it 
also shows that the ~ 7r/5 angle seems to arise in con- 
junction with the Coulomb potential. As we have seen, 
the exponential screening term actually softens this ef- 
fect through the reduction of the pitch angle, but it does 
not suppress it entirely (for small values of n) . From this 
point of view, the chiral symmetry we observe originates 
only from the interplay between the 1/D behavior of the 
Coulomb potential and the cylindrical geometry. Such 
an effect is independent from the specific pattern a and 
therefore we believe it might emerge also in more general 
electrostatic systems over a cylinder. 

The symmetry group of a decorated cylinder falls into 
one of 9 classes, each with a distinct closed subgroup of 
the full cylinder group SO (2) x M The four regions 
in our diagram correspond to the four classes character- 
ized by one-dimensional closed subgroups (the rest being 
zero- or two-dimensional). Namely, up to conjugacy and 
scaling, they are [2^: pure rotations SO(2) x I (fully seg- 
regated phase), continuous rotations and discrete trans- 
lations SO(2) X Z (phase I), discrete rotations and con- 
tinuous translations x R (phase III), and corkscrew 
symmetries Zj, x L (phase II). Z^ is the subgroup of ro- 
tations of the cylinder through angles that are multiples 
of 27r/fc, and L = {(t, i) e SO(2) x M : t e M}. Therefore 
the phase transitions we observe correspond to breaking 
qualitatively different symmetry generators. 

We expect these phases to be relevant for systems at 
the nanoscale, as it can be easily seen by an order-of- 
magnitude analysis of Lq. It is through these length 
scales we believe the functionality of the patterns might 
manifest, such as the induced chirality in our system. 
The advantage of the nanoscale, besides possible practi- 
cal applications, is that electrostatic interactions are still 
important under normal dielectric conditions. Even in 
aqueous systems where charge screening occurs, Coulomb 



interactions still play a vital role through counterions 
and metallic-ion coordination [sO]. Finally, we do not 
explore in this Letter the possibility of a critical point 
where phase I, II or the fully segregated phase meet, 
and we postpone it to future studies. 
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